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ABSTRACT 

We have obtained high resolution spectrograms of smaU scale magnetic structures with the Swedish 
1-m Solar Telescope. We present Doppler measurements at 0'.'2 spatial resolution of bright points, 
ribbons and flowers and their immediate surroundings, in the C I 5380.3 A line (formed in the deep 
photosphere) and the two Fe I lines at 5379.6 A and 5386.3 A. The velocity inside the flowers and 
ribbons are measured to be almost zero, while we observe downflows at the edges. These downflows 
are increasing with decreasing height. We also analyze realistic magneto-convective simulations to 
obtain a better understanding of the interpretation of the observed signal. We calculate how the 
Doppler signal depends on the velocity field in various structures. Both the smearing effect of the 
non-negligible width of this velocity response function along the line of sight and of the smearing 
from the telescope and atmospheric point spread function are discussed. These studies lead us to the 
conclusion that the velocity inside the magnetic elements are really upflow of the order 1-2 kms~^ 
while the downflows at the edges really are much stronger than observed, of the order 1.5-3.3 kms^^. 
Subject headings: Sun: photosphere — techniques: spectroscopic — Sun: atmospheric motions 



1. INTRODUCTION 

In the Sun's photosphere one can observe magnetic 
structures in a range of different scales. From sunspots 
of tens of Mm size to small scale magnetic elements of 
about 200 km or less. Observationally these small scale 
magnetic elements become visible as features brighter 
than the surroundings (leading to the often used term 
"bright points") when the photosph ere is imag ed at suf- 
ficiently high angular resolutions ( Dunn fc Zirker 1973}). 
iTitle fc Bergeil (|1996[ ) showed that bright points can not 
be detected with spatial resolutions worse than about 
0'.'4 due to a smearing out of the contrast because bright 
points are surrounded by darker areas. As observational 
methods have be come more sophisticated, such as adap- 
tive optics (e.g. iRimmeld 120001 : IScharmer et"all l2003bf) 
and post-processing m ethods (e.g. Ivon der Luehelll993t 
Ivan Noort et al.ir2005f) . the observations o f such struc- 
tures have become almost routine. Recentlv lBerger et al.l 
(|2004|) concluded that at 100 km resolution the mag- 
netic elements do not always resolve into discrete bright 
elements. Instead, in plage regions of stronger aver- 
age magnetic field, the small scale magnetic elements 
are found to be concentrated into elongated "ribbons" 
and round "flowers" with a darker core and bright edges. 
These elements are constantly evolving; the merging and 



splitting of "flux sheets" and the transitions between 
the ribbons, flowers, and ni i cro-po res are described by 
iRouppe van der Voort eFall (|2005[ ). The understand- 
ing of these small scale elements is of great importance 
and these structures have therefore be en subjected to 
inten s e research in recent years (e.g., iKiselman et al.l 
20011: iRutten et all [200H I Sanchez Almeida et all |2001| 



Steiner et al. I I2OO1I: iKdlir et al. I 120041: iCarlsson^ et al.l 



2004r ). This has resulted in good theoretical models and 
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a good understanding of why the magnetic elements look 
the way they do in observations. One main conclusion 
of the above papers is that the dominant reason for the 
brightness of bright points is the lower density inside the 
bright point due to magnetic pressure, hence allowing the 
observer to look into deeper layers. These layers are hot- 
ter than the higher layers seen outside the magnetic ele- 
ments because of the temperature increasing with depth. 
At equal geometric height, the temperature inside the 
magnetic element is lower than in the surroundings due 
to the suppression of convective energy transport by the 
magnetic field. This temperature contrast may be large 
enough in larger magnetic concentrations that it offsets 
the effect of the difference in formation height and the 
magnetic element no longer looks bright. Smaller mag- 
netic elements (bright points) and the edges of the larger 
magnetic elements (ribbons, fiowers, pores) are heated 
radiatively by the surrounding hotter, non-magnetic, at- 
mosphere and we get a bright point or a bright edge. 
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Bright points have often been observed in the G-band 
(spectral domain dominated by CH-molecular Hnes at 
4300 A). This is because bright points have higher con- 
trast in this domain due to an increased height difference 
of the optical depth unity height because of the destruc- 
ti on of CH-molecu les in low-density magnetic elements. 

iRimmel^ ()2004l ) presents measurements of the line of 
sight (LOS) velocity field in and around magnetic flux 
concentrations based on narrow-band filtergrams. He 
finds downdrafts at the edge of flux concentrations of 
"a few hundred meters per second." The size of these 
downflow areas is approximately 0'.'2, which indicates 
that they are smaller since this is the approximate spa- 
tial resolution of the observations. Furthermore he con- 
cludes that the downdraft gets narrower and stronger in 
deeper layers and the plasr na within the flux concentra- 
tion is more or less at rest. iRouppe van der Voort et all 
(|2005f) found that magnetic structures such as flowers, 
ribbons and flux sheets have a weak upflow ins ide and a 
sharp downdraft at the edge in accordance with 'R immeld 
(|2004f ). The LOS flow inside flowers and ribbons was 
measured to be about 0-150 ms~^ upflow and the down- 
drafts at the edge to a bout 360 ms~^. 

Bot h lRimmeld ()2004[ ) and lRouppe van der Voort et al.l 
(Hods') use filtergrams at only two spectral positions to 
derive the velocities — adding more spectral sampling 
points would increase the accuracy of the measurements 
and would also make it possible to derive the velocity as a 
function of height. Spec troscopic measurement s of bright 
points were reported bv lLanghans et al] ()2002f ) and they 
derive velocity maps showing downflows clos e to bright 
points. Furthermore iLanghans et al.l (|2004D use spec- 
troscopic diagnostics to measure velocities inside bright 
points, at resolutions of about 0'.'3 they find both upflows 
and downflows inside bright points. 

In this paper, we present spectroscopic data of mag- 
netic elements such as ribbons, flowers, and bright points 
with excellent spatial resolution. Intensity cuts of the 
smallest features and analysis of the spatial power spec- 
trum show that the spatial resolution is better than 0'.'2 
for the spectra used in this work. These spectra are used 
to obtain the LOS velocities in and close to small scale 
magnetic elements. Furthermore we analyze a three di- 
mensional simulation of magneto-convection by solving 
the radiative transfer through the simulation cube in the 
same spectral domain as in the observations. These simu- 
lations are used to understand the effects of the smearing 
along the line of sight (the velocity response function) 
and spatial smearing (the instrument and atmospheric 
point spread function(PSF)) have on the determinations 
of velocities. 

In §[5] we present and discuss the observational program 
and the instrumentation. The data reduction methods 
are discussed in § [3l In § |4] we present the results of the 
Doppler shift measurements in the three spectral lines. 
The analysis of the Magneto-convective simulation, the 
determination of the velocity response functions in the 
diagnostics used, and the discussion of the smearing ef- 
fects this kind of velocity measurements are hampered 
with are presented in § [5] Finally, we summarize our 
results in § [6l 

2. OBSERVING PROGRAM AND INSTRUMENTATION 



The main aim of these observations was to obtain good 
spectra in and close to small scale magnetic structures. 
On 2005 May 13 we observed the NOAA active region 
10759 (7V8.4°, £;i0.5°, /i = 0.97) using a small pore as 
adaptive optics lock point, sec Fig[l]for an example of a 
slit-jaw image. 

The spectra wer e obtained using the&wedish 1-m Solar 
Telescope (SST) (jScharmer e t al.' '20 03aD to gether with 
the adaptive optics system (Scharme r et al.l [2003b) and 
the TRLPort Polarimetric Echelle-Littrow (TRIPPEL) 
spectrograph. As indicated by the name, the TRIP- 
PEL spectrograph is a Littrow spectrograph with an 
Echelle grating using 79 grooves mm~^. The blaze angle 
is 63.43°, the focal length is 1500 mm, and the slit width 
is 25 /im. Simultaneous spectra of two wavelength in- 
tervals, 4566-4576 A (hereafter called the 4571 interval) 
and 5376.5-5388 A (the 5380 interval), were obtained. 
In the present paper we focus on the 5380 interval since 
these spectra were obtained with the shortest exposure 
time and thus the Doppler shifts in the spectral lines 
will be less affected by seeing. The spectrograph should 
always be operated with angles close to the blaze an- 
gle. We observe at 62.73° and in order 47 and 42 for 
the 4571 and 5380 intervals respectively. The theoreti- 
cal slit-limited bandpass for the 5380 interval is W\ — 23 
mA or in terms of effective theoretical spectral resolution 
R « 234000. In all our comparisons with FTS spectra 
and simulations we convolve these spectra to this resolu- 
tion. 

Two Megaplus 1.6 cameras (KAF-0401E Blue Plus and 
KAF-0401 Image sensor) were used at exposure times of 
200 ms and 80 ms for the 4571 interval and 5380 in- 
terval respectively. The two cameras have a quantum 
efficiency of about 35% and 25% at the operating wave- 
lengths. In addition to the spectral cameras we used two 
cameras to obtain slit-jaw images. One of the slit-jaw 
cameras using a 10.3 A bandpass interference filter cen- 
tered at 4572.6 A was slaved to the spectral camera at 
the same wavelength. The other slit-jaw camera, using 
a 9.2 A bandpass interference filter centered at 5321 A, 
had exposure times of 8 ms and was not slaved. In the 
spatial domain the sampling is 0'.'0411 pixel^^. With this 
setup we observed several series of good to excellent see- 
ing with a fixed slit-position. One of these series, with 
a duration of about 17 minutes (10:59:40-11:17:06 UT), 
is excellent both with respect to seeing conditions and 
structures crossed by the slit, see FigdJ In the follow- 
ing all the data presented come from this series and we 
consider data of excellent quality spread over the full 17 
minutes. 

3. DATA REDUCTION 

3.1. Flat fields and dark currents 

The spectrograms are corrected for dark currents and 
flat fields. Flat fields are constructed from 5 series of 50 
images of random scans of the quiet solar disk center. 
We then construct a mean spectrum. To obtain a mean 
spectrum we need to correct for the spectrograph's main 
distortions, namely smile (the curvature of the spectral 
lines) and keystone (different spectral dispersion at dif- 
ferent slit positions). The distortions are determined by 
fitting a 4th-order polynomial to the central part of the 
spectral line and in this way we measure the line core 



3 



3;: h 



o 

in 
o 




1 


Fe I 

1 ^ 










; \ 


r 1 










I 




Fe I 








u 

1 

h 


















■ 




















2 4 6 8 1012 
arcsecs 



5378 



5380 



5382 

A[A] 



5384 



5386 



5388 



Fig. 1. — Slit-jaw image (ie/t) and the corresponding spectrogram (right). The slit-jaw image has been MFBD processed to facilitate the 
identification of small scale structures. The box in the slit-jaw image marks the area which has been magnified in Fig|3] 



position over the slit. This is done for five hnes dis- 
tributed over the spectrogram. The smile is of order 4 
pixels while the keystone is of order 0.15 pixels, thus 
the smile is the most profound aberration. The final flat 
fields are obtained by dividing the total of the 250 images 
by a distortion corrected mean spectrogram. In this way 
we remove both spatial and spectral information without 
losing the fixed CCD dependent signal and without any 
interpolation of pixels in the flat field. 

3.2. Wavelength calibration 

We make a mean solar spectrogram by adding to- 
gether all 250 flat field spectra. The mean aber- 
ration correc ted spectrum is then compared to the 
FTS atlas of iBrault fc Neckell (|1987D . This atlas has 
proved to be well calibrated in wavelength and it 
shows no systematic offset in line shift s with wavelength 
(jAUende Prieto fc Garcia Loped fl998f l. Since the solar 
atlas is corrected for the Earth's rotation, the Earth's 
orbital motion and the Sun's rotation we automatically 
get a corrected spectrogram when we calibrate using the 
FTS atlas. The upper panel of Figl^] gives a good im- 
pression of the accuracy of this calibration. Note that 
in this way our wavelength scale is the same as that of 
the FTS atlas; there is thus no correction for gravita- 
tional redshift. The slit covers a region of almost 38" 
or approximately 28 Mm on the Sun. During the obser- 
vations there is a temporal change of the spectral posi- 
tion on the cameras. This wavelength shift is the same 
in both spectral regions. Part of the shift is periodic, 



with an amplitude of 50-150 ms^^ and a period of five 
minutes. This variation is caused by global oscillations 
averaged over the slit length. In addition there are slow 
trends, probably caused by temperature variations in the 
spectrograph. The frame to frame standard deviation is 
8 ms^^. We correct for this temporal change by ap- 
plying a shift to the wavelength calibration determined 
form the flat field spectra. Since the slit is covering an 
active region, the mean spectrum of each slit is not di- 
rectly comparable to the solar atlas ([Br andt fc Solankil 
Il990l and references therein). iBrandt fc Solankl (J99J3) 
measure the mean convective blueshift in 19 Fe I lines, 
including the Fe I 5379 A line, in quiet Sun and mag- 
netically active regions. They find that the convective 
blueshifts are decreased from 350 ms""'^ in quiet Sun 
to 120 ms~^ in active areas. We determine the tem- 
poral shifts by comparing the line center of the mean 
Fe I 5379 A line in each spectrogram (averaged over an 
active region) with the line center position of the quiet 
Sun atlas. We thus need to compensate for the difference 
in line center convective blueshift between active regions 
an d the quiet Sun. The di fference of 230 ms"^ found 
by iBrandt fc Solankil ()1990r ) refers to an intensity level 
of 60% and not line center. From our observations we 
find that the difference is 30 ms^^ smaller at line cen- 
ter. We thus apply a redshift of 200 ms^^ to the atlas 
before comparing with our mean Fe I 5379 A line when 
determining the temporal shift. 

Instead of correcting each spectrum for smile and key- 
stone we use the determined aberrations to determine 
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Fig. 2.— Upper: FTS atlas of the Fe I 5379.6 A line (solid) 

and the mean spectrum (dashed- dotted) . The scattered 
light corrected mean spectrum is showed with triangles. 
Lower: Intensity level (intensity relative to the mean con- 
tinuum intensity) of our quiet Sun spectrum (mean flat 
field spectrum) versus the intensity level of the FTS at- 
las (crosses) and a least-square linear fit to these data 
(solid). 

a separate wavelength scale for each slit position. This 
procedure avoids interpolation of the spectrogram data. 

3.3. Scattered light 

When the mean quiet Sun spectrum is plotted against 
the FTS atlas, see Figl^] upper panel, it is clear that the 
mean spectrum has shallower lines than the FTS atlas. 
This discrepancy is due to scattered light, mainly from 
the diffuse scattering in the spectrograph. In the lower 
panel in Figl^] we have plotted the intensity level of the 
mean spectrum against the FTS intensity level. We use 
a least-square linear fit to these data points. Assuming 
a constant level of diffuse scattered light, the offset and 
the new continuum calibration of the mean spectrum is 



given by this fit. The amount of constant scattered light 
is 6 %. This level is subtracted from the spectra and 
the resulting spectra are then divided by a scaling factor 
found from the slope of the curve in the lower panel of 
FigO Again we refer to Figl2] upper panel to see the 
effect of this co rrection of the s catter ed light. 

According to ISobotka et all (|1994D . "the mean inten- 
sity in abnormal granulation is equal to that of quiet 
granulation," so we expect the mean spectrum from the 
fiat fields to have the same mean count as we have in 
each of the mean spectrograms when we do not include 
the pores. This turns out not to be the case. There are 
many reasons for this discrepancy, different atmospheric 
conditions, the changing path length of the light ray in 
the atmosphere due to the time of the day and changing 
intensity on the Sun due to solar oscillations (p- modes). 
All these effects are corrected for by introducing a scal- 
8 ing factor which gives the same mean solar continuum 
intensity in the mean spectrum of each spectrogram as 
in the mean solar atlas. We multiply with this scaling 
factor before we correct for scattered light. 

3.4. Slit-jaw images 

The slit-jaw images are corrected for dark currents and 
fiat fields and post-processed usi ng the mu lti frame blind 
deconvolution (MFBD) method (Ivan Noort et al. 2005). 
The intensity on the slit in the slit-jaw image is corre- 
lated to the mean intensity in the corresponding spec- 
trogram to obtain the relative scaling and offset between 
the spectrogram and the slit-jaw image. The aligned slit- 
jaw images are used to identify the structures crossed 
by the slit, see Fig[TJ Note that the slit-jaw image has 
shorter exposure time than the spectrum and has been 
post-processed — it is thus sharper and has a more nar- 
row point-spread-function than the spectrogram. 

4. OBSERVED DOPPLER SHIFTS 

The 5380 interval contains the deep forming 
(|Livingston et all I1977D C I 5380.3 A line, see FigdJ 
Furthermore it contains a well suited quite strong Fe I 
line at 5379.6 A and a weaker Fe I line at 5386.3 A. 
These two iron lines have been singled out by other au- 
thors as good both in the sense that they are blend free 
and thus are well suited to measure the line asymme- 
try (jPravins et al.l 1 19811] and that their central wave- 
lengt hs are well kn o wn ([Allende Prieto &: Garcia Loped 
Il998f) . iNave et al.l ()1994 f) give the central wavelength 
of these Fe I lines, see Table [TJ and they claim the er- 
ror in the central wavelength of these two lines to be 
less than 1.25 mA and 2.5 mA respectively, which cor- 
responds to uncertainties in velocities of 70 ms~^ and 
140 ms ~^. The central wavelength of the C I line is 
given bv I Johansson! (|1966D and he quotes the uncertainty 
to be less than 0.02 A, which corresponds to 1115 ms^^. 
This high uncertainty makes it difficult to use this line 
to determine velocities on an absolute sc ale. The Fe I 
line a t 5383.3 A was also picked out by iDravins et al.l 
(|T98lh . but we do not use this line due to a blend in the 
blue wing. We use the laboratory wavelengths as ref- 
erence wavelengths for the two Fe I lines, while we use 
a reference wavelength which gives us a more realistic 
convective blues hift for the C I line. Using the wave- 
length given bv iJohansso^ (fl966il . 5380.337 A, we get 
a convective blueshift of 1465 ms~^. Instead we use a 



5 



reference wavelength which is in agreement with a 3D 
magneto-convective simulation, see § [S] The simulation 
gives a convective blueshift of 664 ms~^, which indi- 
cate a convective blueshift in non-magnetic regions of 
864 ms~^, which corresponds to a reference wavelength 
of 5380.3262 A. This wavelength is well within the given 
uncertainty of 0.02 A. The line-core wavelengths in the 
FTS atlas, the reference wavelengths, used to calculate 
the LOS velocities, and the corresponding blueshifts are 
shown in Table [TJ Notice that we have corrected for a 
gravitational redshift of 633 ms~^. The Lande g-factors 
of the Fe lines are 1.10 and 1.17, for the C-line it is un- 
known but greater than zero. 



Central wavelengths and convective blueshifts 



Atomic 


Atlas 


Reference 


Velocity 


species 


[A] 


[A] 


[ms^-"-] 


Fe 1 


5379.5800 


5379.5740 


-307 


C I 


5380.3221 


5380.3262 


-864 


Fe I 


5386.3349 


5386.3341 


-589 



TABLE 1 



Since the C I line is very shallow we only determine the 
total line shift for this line, while the Fe I lines are well 
suited for determining the Doppler shift at different in- 
tensity levels by determining the bisector of the line. To 
facilitate a comparison of velocities in different features 
we use the same absolute intensity levels in the various 
features (we thus do not use the local continuum of that 
feature). The intensity levels used are relative to the av- 
erage continuum level of the aver age Sun at disk center . 
This continuum value, as given bv lNeckel fc L^(fl98l . 
and the values at other intensity levels are given in Table 
m both as absolute intensities and as radiation tempera- 
tures. The latter is indicative of the gas temperature at 
the level of formation of the radiation at this intensity 
level. 



Intensity level, specific intensity, 



and radiation temperature. 


IL 


I. 


RT 




[ergcm~^str~^ s~^ Hz"-"^] 


[K] 


1.00 


3.70E-05 


6293 


0.95 


3.51E-05 


6219 


0.90 


3.33E-05 


6143 


0.85 


3.14E-05 


6064 


0.80 


2.96E-05 


5983 


0.75 


2.77E-05 


5899 


0.70 


2.59E-05 


5811 


0.65 


2.40E-05 


5720 


0.60 


2.22E-05 


5625 


0.55 


2.03E-05 


5524 


0.50 


1.85E-05 


5418 



TABLE 2 



The total line shift in the C I line is calculated using a 
Gaussian fit to the line. When we calculate the velocities 
we correct for gravitational redshift. 

Our slit crosses from two to three ribbons, depending 
both on the temporal evolution of the photosphere as well 
as the changes in slit position due to differential seeing. 
It must be emphasized that we do not have magnetic in- 
formation so all classification has to be done purely by 
visual inspection. The kind of ribbon-structures we see 



has been shown to be magnetic (|Berger et al.ll2003 ). A 
flower/ribbon and a typical granule can be seen in FiglHl 
This Figure also shows velocities and the continuum in- 
tensity at the slit position. It is clear that we have down- 
flows at the edge of the flower/ribbon (the two arrows), 
which are increasing with increasing intensity level in 
the line (formation deeper in the atmosphere). Further- 
more, the velocity in the granule is somewhat high (the 
mean measured value in our observations is about 1.2- 
1.7 kms""'^), with about 1.4-2.0 kms~^ upflow. Finally, 
we see that in the center of the flower/ribbon we measure 
the velocities to be more or less zero. Intensity minima 
between magnetic structures and granules are found by 
visual inspection, these are hereafter called "intergranu- 
lar lanes close to magnetic structures" (IGM). The peak 
velocity close to the intensity minimum in the IGM is 
measured, this peak can vary slightly from the inten- 
sity minimum, but not more than 0'.'15, see Fig|3l This 
peak velocity is measured for up to 51 different IGMs, 
the velocities and the standard deviations are shown in 
Table [31 The pixel-to-pixel standard deviation in the 
velocity determination is 40 ms^^. 



Velocities in IGMs 



Atomic 
species 


IL 


Velocity 

[ms-i] 


[ms-i] 


Samples 


C 1 5380 




191 


227 


51 


Fe I 5379 


0.90 


755 


326 


51 




0.85 


563 


225 


51 




0.80 


459 


188 


51 




0.75 


407 


179 


51 




0.70 


364 


173 


51 




0.65 


334 


183 


51 




0.60 


314 


182 


47 




0.55 


283 


179 


35 


Fe I 5386 


0.90 


410 


273 


51 




0.85 


278 


193 


51 




0.80 


269 


216 


44 



TABLE 3 



arcsecs 

2 4 6 8 



Relative intensity 




-2000 -1000 1000 2000 

Velocity [ms~'] 

Fig. 3. — One of the snapshots shows the spectrograph's slit 
crossing a flower/ribbon (20'.'8-22^'7) and a granule {18'.'5-20('0). 
The velocities are shown as measured in the Fe I 5379.6 A line at 
intensity level 0.85 (red) and 0.65 (green), in the C I line [blue) and 
in the Fe I 5386.3 A line at intensity level 0.85 (yellow). Positive 
velocities mean downflow. The black line shows the relative con- 
tinuum intensity in the spectrogram. The two IGMs are marked 
with arrows. 
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A statistical analysis of the minimum LOS velocity 
within flowers and ribbons (Ribbon centers or RC) gives 
us velocities around zero ms~^ in the Fe I 5379 A 
line, while the the C I line shows an up flow of about 
300 ms~^. The Fe I 5386 A shows a small upflow inside 
RCs with a mean of about 200 ms~^. The standard de- 
viations are of the order 100-200 ms^^ in the Fe I lines, 
while the C I line has a standard deviation of about 
350 ms~^. These results are hampered with relatively 
high uncertainties of a few hundred ms~^. We therefore 
want to investigate the effects of smearing along the line 
of sight and spatially on these velocity measurements. 

We have also measured velocities inside bright points 
but since the spatial extent of each bright point is small 
they are difhcult to identify. Nevertheless we visually 
classify 15 bright points and measure the velocities at 
intensity maximum. We measure a mean downflow of a 
few hundred ms~^, but the velocities inside bright points 
vary a lot, and we also measure up flows inside some of 
the bright points. The standard deviation is quite high, 
about 0.5 kms^^. Due to the small spatial extent of 
each bright point, the measured velocities are strongly 
affected by the surrounding intergranular velocity. 

The effect of a non-zero g-factor is increased broaden- 
ing in magnetic areas. For a non-symmetric profile this 
broadening may also introduce a change in the deter- 
mined bisector. We have quantified this effect by apply- 
ing a Zeeman splitting corresponding to a magnetic field 
strength of 1000 Gauss to the observed profiles and rede- 
termined the bisector. The change in measured velocity 
is less than 40 ms~^. 

5. SIMULATIONS 

Since the intensity level can be considered as a height 
parameter it is in principle possible to obtain the ve- 
locity as a function of absolute physical height. In or- 
der to explore this possibility we need a good model 
with which we can calibrate the velocities with height. 
One very suitabl e model is the magnet o-convective sim- 
ulation used by ICarlsson et al.l (|2004[ ). In this realis- 
tic simulation we can see bright points as well as more 
elongated magnetic structures. To tailor this simulation 
to our observational setup we use 12 atmospheric snap- 
shots, includiri g the same snapshot as the one used by 
ICarlsson et al.l (|2004l ). with one minute time difference 
between the snapshots. We solve the equations of radia- 
tive transfer in LTE for these 3D atmospheres, using the 
MULTI code (,Carlsson.4986, ). for 1364 frequency points 
spanning from 5377.6 A to 5389.4 A. A Hue list of 249 
lines is used. The line data was retrieved from the Vi- 
enna Ato mic Line Da t abase (V ALD) (Kupka et al. 1999; 
iPiskunov et al.l Il995l : iRvabchikova et al.l Il999|[ ). with a 
few exceptions; the C I oscillator strength is taken from 
iHibbert et all (|1993D and an Fe line at 5382.5 A was re- 
moved due to obvious errors in its central wavelength and 
oscillator strength. 

A temporally and spatially averaged mean spectrum 
over 12 minutes and 8'.'27 x 8'.'27, totally over 750000 
spectra, is constructed. Since the simulation is cover- 
ing an active region we expect the mean spectrum to 
have less blueshift, have shallower lines and to have 
more vertical lower bisectors (intensity level less than 0.7) 
than s imilar mean spectra obta ined in non-magnetic re- 
gions (|Brandt fc Solankil Il990l . and references therein). 
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Fig. 4. — From the FTS solar atlas, the total magnetic simula- 
tion and the total non-magnetic simulation we derive the bisectors 
in the two Fe I lines. The mean bisectors in the Fe I 5379 line 
{blue) and the Fe I 5386 A line (red) in the FTS atlas (solid), the 
non-magnetic simulation (dashed), and in the magnetic simulation 
(dashed- dotted) are shown. The shift in the Fe I 5386 A line indi- 
cate that there might be an error in the laboratory wavelength of 
about 100 ms~^. Furthermore the laboratory wavelength of the 
Fe I 5379 A seems to be well determined. Also notice the reduced 
blueshifts in the magnetic simulation of approximately 200 ms~^, 
which are in good agreement with observational results. 

For reference we therefore solve the equations of ra- 
diative transfer in LTE in three non-magnetic snap- 
shots with 1 minute cadence, totally over 190000 spec- 
tra. Even though the temporal sampling is quite low 
we believe that the difference from a bigger temporal 
sample will be within 1 mA or approximately 50 ms~^. 
In the non-magnetic simulation we have deeper strong 
Fe lines than the atlas, this is d ue to NLTE effects 
(jShchukina fc Truiillo Buenoll2001h . The Fe hues in the 
mean magnetic spectrum do not show any reduced line 
depth. This might be because the magnetic field in the 
simulation is too weak but also because of NLTE ef- 
fects. The convective blueshifts are in general w ell re- 
produ ced in the non-magnetic simulation ( Asplund et al.l 
l2000f ). One difference is that the upper part of the bisec- 
tor shows larger blueshift than the atlas, see FiglH We 
therefore use the blueshift in the magnetic simulation as 
reference wavelength for the C I line since this simula- 
tions seems to give a more realistic convective blueshift 
in this line, see § |4l The magnetic simulation shows a 
blueshift that is about 200 ms^^ smaller than the non- 
magnetic simulat ion, very similar to what was found ob- 
servationally by (^ Brandt fc Solankil |1990[ ). Furthermore 
the lower part of the bisector in the magnetic simulation 
is more vertical than in the non-magnetic simulation, also 
similar to observations. These are indications that the 
simulations are reproducing the velocities in a realistic 
manner. 

When we measure the velocity in a spectral line at a 
given intensity level there are two smearing effects that 
will influence the measurement. The first effect is an in- 
trinsic property of the line formation mechanism, namely 
the finite width of the response function. The response 
to a given velocity at a given height in a spectral line 
at a given intensity level is given by the respons e func- 
tion to velocity at this intensity level (iMagainI |1986X 
We calculate the response functions numerically using 



a step function in the vel ocity, very similar to the way 
iFossum &: CarlssonI (|2005[ ) calculate their response func- 
tions to perturbations in temperature. Typical response 
functions will have a pronounced peak but with a sig- 
nificant width of a few hundred kilometers, see FiglSHH] 
for some examples. In Table |4] we give average response 
heights (given as the first order moment) in six different 
solar features (see Fig l5] and Figl7]for definition). These 
response heights are calculated in the line core (LC) and 
in the close to continuum line wing (CLW) for the two 
Fe I lines while the moment of the total line (TL) was 
used for the C I line. 

First order moment of the response functions. 
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TABLE 4 



The width of the response functions means that we do 
not measure the real plasma velocity but rather a height 
smeared velocity. The difference between the response 
function weighted velocity and the velocity at monochro- 
matic optical depth unity can be several hundred ms~^, 
see FiglT] It is thus important to be aware of this effect 
when one talks about Doppler shifts measured in spec- 
tral lines. The response function to velocity in the IGIVI 
is quite narrow at all intensity levels, while the response 
function in the RC is much wider, see FiglS] The ef- 
fect of the response function smearing is, however, much 
higher in the IGIvI since the actual atmospheric velocity 
is changing much more quickly with height in IGlVIs than 
in RCs. This explains why we measure almost the exact 
plasma velocity in the RC in FigiT] while we measure 
lower velocities in the IGIVI. 

The second effect is the spatial smearing that origi- 
nates from atmospheric seeing, the diffraction pattern of 
the telescope and scattering in the atmosphere, the tele- 
scope and the spectrograph. Even if we could determine 
the total point-spread function (PSF) of the atmosphere- 
telescope-spectrograph system, we can't correct for these 
effects in the measured spectra since we have spectral in- 
formation only along the one dimension of the slit. An- 
other approach is to model the PSF and convolve the sim- 
ulations with such a model to quantify the effects of the 
PSF on the measurements. The PSF of a partially AO 
corrected image for short exposures is usually assumed 
to consist of a diffraction core and a seeing halo from un- 
corrected or partially corrected modes. We assume that 
we can use this model for our spectra even though 80 ms 
is longer than what is normally accepted as short expo- 
sure. Furthermore we assume that the seeing halo has the 
form of a Lorentzian ()Nordlundlll984f ). Our model PSF 
thus has two free parameters: the width of the Lorenzian 
and the fraction of the Lorentzian component compared 
with the Airy core. A more detailed model of the PSF 
chould not be obtained since we lack the statistics neces- 
sary for the determination of more than two parameters. 
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Fig. 5. — Some typical response functions to velocity. The re- 
sponse function of the total C I line Doppler shift (red), and of 
the line core Doppler shift for the Fe I 5379.6 A line (blue) and 

Fe I 5386.3 A line {green) are shown. The different 
panels show from upper left to lower right: Ribbon cen- 
ter (RC), Intergranular lane close to IMagnetic element 
(IGM), Bright point (BP), Granule (G), Large IVIagnetic 
Concentration (LIVIC), Inter granular lane (IG). The dif- 
ferent points are defined in Figl?! the BP , G, LIVIC, and 
IG are the same points as those used bv ICarlsson et all 
(|2004). 

The observable is the continuum intensity distribution 
as measured in the spectra, excluding the pore. We as- 
sume that the intensity distribution in the simulations is 
realistic and compare the normalized distribution func- 
tion of the intensities in the smeared simulation with the 
observed distribution function to determine the parame- 
ters. We find a good fit with a Lorentzian with FWHIM 
of O'.'T having a fraction of 0.1 at the center of the PSF, 
see FigH The Strehl ratio of this PSF is 0.15 which is 
low compared with what one would expect from the AO 
system of the SST. One should remember here that we 
have no separate component to describe the scattering of 
the Earth's atmosphere, which is known to be very wide, 
and the Strehl ratio obtained includes this effect. It is 
also important to point out that a range of parameters 
give acceptable fits to the observed intensity distribu- 
tion function — wider Lorentzians with smaller fraction 
gives almost the same effect on the normalized distribu- 
tion functions as a narrower Lorentzian with larger frac- 
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Fig. 6. — Height variation of the response functions in the Fe I 
5379 A Hne in the Ribbon center (RC) (red) and the Intergranular 
lane close to Magnetic element (IGM) (blue). The intensity levels 
are 0.4 (solid), 0.6 (dashed) and 0.8 (dot-dashed). These intensity 
levels roughly correspond to the line core (LC), the mean line wing 
at the mean intensity level between the LC and local continuum, 
and the line wing close to continuum (CLW) respectively. 



granules have large upflows, low intensity intergranular 
lanes downflow. When the PSF gives the observed inten- 
sity distribution function this also implies a similar effect 
on the velocity distribution function. Before comparing 
with observations, the spectral domain is convolved with 
the theoretical resolution of the TRIPPEL spectrograph. 
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Fig. 7. — Granule and a ribbon as they appear in the simulations 
are shown in this simulated slit-jaw image. The plasma velocity at 
the height where we have an atmospheric temperature equal to the 
radiation temperature at an intensity level (IL) of 0.65 is given in 
yellow, the non-convolved measured velocity at IL=0.65 in the Fe I 
5379 A line is given in red, while the green and blue lines show the 
convolved measured velocity at IL=0.85 and IL=0.65 respectively 
in the Fe I 5379 A. The relative slit continuum intensity is shown 
in black. In Fig[5] we show response functions for different solar 
structures, these structures are defined in this figure: -|- shows RC, 
o shows IGM, X shows BP, O shows G, □ shows LMC, and A 
shows IG. 

tion. To check the sensitivity of the results to the PSF, 
we have performed tests both with a wider central core 
than the Airy disk and with different combinations of 
the width and fraction of the Lorentzian. As long as the 
normalized distribution function is close to the observed 
one we do not see significant changes in the resulting 
velocities (in comparison with the standard deviations 
shown in FiglS]). At the same time, the velocities deter- 
mined from the convolved simulations are very different 
from the velocities in the unconvolved simulation. This 
may seem contradictory but is due to the fact that in- 
tensities and velocities are correlated — high intensity 



Fig. 8. — To mimic the atmospheric seeing we convolve the sim- 
ulations with a linear combination of a Lorentzian and the Airy 
disk. The resulting NDFs for the observed spectra (solid), the un- 
convolved simulation (dashed), and the convolved simulation (dot- 
dashed) are shown. 



Central wavelengths 
and convective blueshifts in the simulations 
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1.0 F 



0.9 r 



a; 
> 
CO 

- 0. 



0.7 ^ 



0.6 r 



0.5 t 







It- / 








r '' 77 









200 



400 



600 



800 1000 1200 



Fig. 



Velocity [ms ] 
9. — Measured velocities in fhe IGMs m the observations 



(blue) and in the simulations (red) measured at different intensity 
levels for the Fe I 5379 A line (solid), the Fe I 5386 A line (dashed) 
and the C I line (asterisks). The standard deviations for the Fe I 

5379 A hne are shown with error bars in the respective 
colors. 

We use the laboratory wavelengths as reference wave- 
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lengths for the Fe I hnes, while we use the same wave- 
length as in the observations for the C I line. Further- 
more we calculate the central wavelength in the simu- 
lated atlas and the corresponding blue shift. The results 
are shown in Table [5] (with negative velocity meaning 
blueshift). 

Now we can measure the velocities in the IGMs in the 
simulations. This is done in a similar manner as we mea- 
sured the velocities in the IGMs in the observations, the 
results are plotted together with the observed results in 
FigO The absolute velocities in the simulations and the 
observations differ with 200 ms~^ or less, which is amaz- 
ingly close when one consider the uncertainty in the ob- 
servations of about 100-200 ms~^ and that the standard 
deviation is about 200 ms~^. Note that we believe that 
the Fe I 5386 A hue has an error of about 100 ms ^ in 
the laboratory wavelength. This will increase the mea- 
sured velocities with the same amount. Furthermore we 
see an increase in measured velocity with intensity level 
for both the observations and the simulations. This in- 
crease is somewhat lower in the simulation than in the 
observations. Note that the C I line gives substantially 
lower downflow values than the iron lines at compara- 
ble intensity level. This is true both in the observations 
and in the simulations. In the latter, the determined ve- 
locity is completely independent on the reference wave- 
length used so an error in the reference wavelength is not 
the explanation for the low velocities in the observations. 
The explanation is rather the differences in the response 
functions and the spatial smearing due to the seeing. The 
difference between the response functions in the IGM is 
minimal, but the difference between the response func- 
tions in granules is substantial. In the granule the C I 
line is formed deeper in the atmosphere, and hence this 
line will show a larger upflow. When the spatial smearing 
is applied this difference will lead to a smaller downflow 
in IGMs in the C I line, as observed. It is clear that the 
simulations and the observations are not in disagreement, 
even though there is room for improvements. 

We also measure the lowest velocities within the mag- 
netic elements in the simulation. The velocities measured 
range from about 200 ms~^ upflow to about 200 ms^^ 
down flow. The standard deviation is 250-300 ms 
Due to the lack of good sampling of bright points in the 
simulations we do not include any statistics on bright 
points in the simulations. 

6. CONCLUSIONS 

We have presented spectrograms of small scale mag- 
netic structures with excellent spatial resolution. These 
spectrograms have been used to measure the LOS veloci- 
ties in several different magnetic structures. At the edges 
of the ribbons and flowers we measure downflow with an 
increasing magnitude with depth in the atmosphere. The 
magnitude of the downflow is about 200-700 ms~^ with 
increasing velocity deeper down in the atmosphere. In- 
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side the ribbons and flowers we measure a small upflow. 
The upflows are usually in the range 100-300 ms~^. Our 
investigation of the velocities inside bright points shows 
both up and downflow, nevertheless with a mean down- 
flow of a few hundred m s^-'^ . In most cases the velocity is 
closely related to the intergranular velocity, which indi- 
cate a strong presence of scattered light in the measured 
velocity. 

Furthermore we have investigated the processes that 
are involved in the smearing of the measured velocities. 
The results from this work can be summarized as follows: 

1. We have calculated response functions and response 
heights for several different solar features. The width of 
the response functions is usually a few hundred km. 

2. We have convolved the simulations with a simple 
model of the point spread function of the observations. 
Due to the strong horizontal gradients in the velocity 
fleld in and close to small scale magnetic features, this 
spatial smearing has a strong effect on the measurements: 
the measured velocities in the intergranular lanes close 
to magnetic elements are decreased with about 1 kms~^. 

3. Since the simulations are reproducing the observa- 
tions reasonably well we believe that the plasma veloci- 
ties in the simulations are close to the real solar plasma 
velocities. In particular the velocities in and close to 
larger magnetic structures such as ribbons are well re- 
produced. The velocity in the simulation in intergranular 
lanes close to magnetic structures increases from about 
1.5 kms^^ at 150 km height to 3.3 kms~^ at km. 

The velocity inside bigger magnetic structures such as 
ribbons in the simulations is about 1-2 kms^^ upflow, 
with the velocity changing only a few hundred ms~^ 
throughout the formation height between about —100 km 
and 250 km. 
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